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Abstract 

Performance of quadrature amplitude modulation (QAM) schemes is analyzed with dual-hop decode- 
and-forward (DF) relaying systems over mixed rj-fi and fading channels. Closed-form expressions 
are obtained for the average symbol error rate (ASER) for general order rectangular QAM and cross 
QAM schemes using moment generating function based approach. Derived expressions are in the form 
of Lauricella’s (•), $^"^ (•)) hypergeometric functions which can be numerically evaluated using 
either integral or series representation. The obtained ASER expressions include other mixed fading 
channel cases addressed in the literature as special cases such as mixed Hoyt, and Rice fading, mixed 
Nakagami-m, and Rice fading. We further obtain a simple expression for the asymptotic ASER, which 
is useful to determine a factor governing the system performance at high SNRs, i.e., the diversity 
order. Additionally, we analyze the optimal power allocation, which provides a practical design rule 
to optimally distribute the total transmission power between the source and the relay to minimize the 
ASER. Extensive numerical and computer simulation results are presented that confirm the accuracy of 
presented mathematical analysis. 
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I. Introduction 

Cooperative eommunieation has reeeived eonsiderable attention for several emerging wireless 
network arehitectures, sueh as eellular networks, loeal area networks and heterogeneous networks 
to eounter the effeet of multipath fading whieh degrades the symbol error rate (SER) performanee 
m-m. The eooperative eommunieation offers signifieant inerease in eapaeity and multiplexing 
gain in the aforementioned wireless networks where user systems and nodes in a wireless 
network share their resourees and ereate eollaboration through distributed transmission. The 
motivation behind the eooperative eommunieation is to provide virtual multiple-input multiple- 
output (MIMO) support for a user system whieh eannot aeeommodate multiple antenna due to 
size, eomplexity, power eonsumption and many other eonstraints. In a eooperative eommunieation 
system, user terminals work as information sourees as well as relays. The two main relaying 
methods used in cooperative communication scheme are amplify-and-forward (AF) and decode- 
and-forward (DF). In the AF case, the relay amplifies a signal transmitted by the source and then 
re-transmits the same to the destination. In the DF case, the relay decodes a signal transmitted 
by the source and then encodes and transmits to the destination. 

In recent years, there has been an increased interest on the performance analysis of dual-hop 
cooperative communication networks where links are subject to asymmetric fading conditions 
[|6l - [fTTll in which the propagation in one of the links is dominated by a strong line-of-sight 
(FOS) component relative to other links in the channel IfTSl . The performance analysis of dual¬ 
hop AF relay networks over mixed fading conditions are available in [jH- [|T3l and the analysis 
for dual-hop DF relay networks over mixed fading conditions are studied in [fTSll - [fTTll . In IfTTll . 
authors present exact SER expressions of M-ary phase shift keying (MPSK) scheme for DF 
relay system over n-fi and 77-/2 and mixed and 77-/2 fading channels. However, no closed- 
form SER expression of MPSK in mixed K-jj, and 77-/2 fading channels has been provided. 
There are a number of works reported on the performance of the DF cooperative diversity over 
fading channels focusing on the capacity, diversity gain, and outage behavior of the system. 
Past few years has seen a research interest on the investigation of error performance of various 
modulation schemes for the DF cooperative diversity over general fading channels. In [fT9ll and 
ll20l . SER performance for DF cooperation is examined over Rayleigh fading channels. The 
SER performance of different modulation schemes for dual-hop DF cooperation systems have 
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been derived over Nakagami-m fading channels in IfUl - If24]| . In 11211 . the end-to-end bit error 
rate (BER) performance of cooperative diversity networks using DF relaying over independent 
non-identical flat Nakagami-m fading channels is available. The SER expressions of MPSK 
and quadrature amplitude modulation (QAM) schemes for DF cooperative communications over 
Nakagami-m are derived in ll22l . In ll2^ . average BER performance of DF cooperative systems is 
investigated for binary PSK (BPSK) signals in Nakagami-m fading channels for integer values of 
m. Error probability of opportunistic DF relaying in Nakagami-m fading channels with arbitrary 
m is given in [l24l . The authors in ll25ll present the outage and error rate performance analysis of 
the DF cooperative diversity system with OSTBC over spatially correlated Nakagami-m fading 
channels for integer values of m. The error performance of MPSK and square QAM (SQAM) 
signals for dual-hop DF cooperation systems have been investigated over two-wave with diffuse 
power fading channels in Il26l . In ll27ll . authors analyze the performance of SQAM schemes 
with DF relaying system over Hoyt (Nakagami-g) fading channels. In Il2^ . the exact SER of 
M-ary PSK for multihop communication systems with regenerative relays is provided. In [[29ll . 
the outage probability (OP), SER, level crossing rate and average outage duration are derived 
for multihop DF cooperation systems over Generalized-iT fading channels. The OP and BER of 
Gray-coded rectangular QAM (RQAM) signals for multihop DF cooperative systems have been 
analyzed over rj-jd fading channels in ll^ . A closed-form expression for the OP of dual-hop DF 
relaying in dissimilar Rayleigh fading channels is derived in llSTI . The authors in |[^ and |[3^ 
are investigated the OP for the dual-hop DF protocol over Nakagami-m fading channels. The 
OP and outage capacity of dual-hop DF relaying system over and K-jj, channels are given 
in [f34l and (SSI, respectively. 

It is well known that the classical small-scale fading models (e.g.. Rice, Hoyt, Nakagami-m) 
are not able to accurately characterize the small-scale variation of the fading signals, particularly 
at the tail probability |[^ . Yacoub addressed this problem suggesting two new distributions for 
fading models namely and n-fi which better accommodates practical situations for which 
classical distributions are a poor fit. The 77-/1 distribution accurately models small-scale fading 
for various non-FOS (NFOS) conditions including the Hoyt, and Nakagami-m distributions as 
special cases. On the other hand, the n-jj, distribution is a generalized distribution for modeling 
a great variety of EOS channels and includes Rice and Nakagami-m distributions as special 
cases. Some available works on performance analysis issues of cooperative systems operating in 
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these two generic fading channel models are available in lfT3]| . IfTTlI . OOl . [|34l . [1^ where the 
performance of AF and DF cooperative systems are investigated. 

QAM is a promising modulation scheme for a digital multimedia transmission in wireless 
communications due to its high bandwidth efficiency [ITTlI . RQAM, SQAM and cross QAM 
(XQAM) are well known QAM schemes among which RQAM is a generic modulation scheme 
which includes SQAM, binary PSK (BPSK), orthogonal binary frequency-shift keying, quadra¬ 
ture PSK (QPSK) and multilevel amplitude shift-keying modulation schemes as special cases. 
XQAM is an optimal QAM constellation for odd number of bits per symbol as it has low average 
symbol energy than RQAM [l38ll . XQAM finds application with constellations from 5 bits to 15 
bits with usage in asymmetric digital subscriber lines and very high speed digital subscriber lines. 
Specifically, 32 and 128-XQAMs are employed in digital video broadcasting-cable. Moreover, 
XQAMs have been found to be useful in blind equalization [|^ and adaptive modulation [|4^ . 

m. 

Despite many research works dealing with the DF cooperative diversity during the past few 
years, to the best of our knowledge, closed-form expressions for the average SER (ASER) of 
general order RQAM and XQAM schemes with dual-hop DF relaying systems over mixed rj-fx 
and K-fx fading channels is not available in literature. In this paper, we derive novel closed- 
form expressions for the ASER for general order RQAM and XQAM schemes with dual-hop 
DE relaying systems operating over mixed Tj-fx and K-/i fading channels. The obtained ASER 
expressions contain special functions such as Eauricella’s hypergeometric func¬ 

tions which can be easily evaluated numerically using their finite integral or converging infinite 
series representation [l42l - Pl6l . It is worth mentioning here that the derived expressions include 
the ASER expressions for mixed Hoyt, and Rice fading, mixed Nakagami-m and Rice fading, 
and non-identical Nakagami-m fading as special cases. We further derive a simple expression for 
the asymptotic ASER, which is useful to determine a factor governing the system performance 
at high SNRs, i.e., the diversity order. Additionally, we obtain the optimal power allocation 
factor, which provides a practical design rule to optimally distribute the total transmission power 
between the source and the relay to minimize the ASER. 

The rest of the paper is organized as follows. Section II deals with system and channel 
models and the ASER expressions are derived in Section III. In Section IV, numerical results 
and discussion are given. The paper is concluded in Section V. 
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II. System and Channel Models 

We consider a three-node system model ^22^ . This model finds a typical application of 
cooperative communications for the uplink of cellular wireless systems. This type of cooperative 
communication is also suitable for handsets equipped with single antenna for achieving transmit 
spatial diversity for link-quality improvement. In this system model, the source node s sends 
its information to the destination node d on two consecutive time slots. In the first time slot, 
node s broadcasts its symbol to node d and relay node r. In the second time slot, only node 
r, if decodes successfully, forwards the received symbol to node d. In this paper, we assume 
that the relay can successfully decode if the SNR level of the sr link is above a predetermined 
threshold. In a single relay dual-hop DF relaying system, node i sends its information to node 
j over ij link, where i G {s,r}, j G {r, d} and ij G {sd, sr,rd}. When unit energy symbol x 
is transmitted from node i, the baseband signal received at node j can be given as 


Vij 




( 1 ) 


where P, is the power of the transmitted signal at node i, aij denotes fading channel coefficient 
of ij link which is modeled as either //-/i or K-/r distribution, is additive white Gaussian 
noise (AWGN) with Nq variance at node j. At node s the symbol is transmitted with power 
Pi = Ps. If node r successfully decodes the received symbol, the symbol is transmitted with 
power Pi = Pr, otherwise node r remains idle, i.e. P^ = 0. The total power P of this dual-hop 
DF relaying system is given by P = P^ -f P^, 0 < P^, Pr < P. Finally, node d combines the 
signals of both the time slots according to maximal ratio combining (MRC) technique. 

The instantaneous SNR, 7 ^ of ij link is defined as 7 *^ = a^jPi/No. If a link is subjected to 
rj-jd fading, the probability density function (PDF) of 'jij can be expressed as Il36l 


Ui ( 7 ) 


- - - J - — exp 


2/i j j hij 7 

lij 


IJ-ij- 


‘2‘^ijHij'y 

lij 


( 2 ) 


where r(-) is the gamma function, I^{-) is the modified Bessel function of the first kind and 
uih. order, Hij denotes the number of multipath clusters, and both hij and Pdij are functions of 
rjij. In dH), 'jij = E[ 7 jj] = ilijPi/No denotes the average SNR of the ij link , where E[-] is 
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the expectation operator and VLij = Eiafj] is the variance of aij. The distribution in @ has 
been described for two types of physical models (Formats) depending on the way the parameter 
r]ij is defined. In this paper, we consider only Format 1. In Format I, 0 < rjij < oo is the 
ratio of the power of inphase and quadrature phase components of the scatter wave signals in 
each multipath cluster, with hij = (2 + + riij)/A and Hy = — ? 7 jj)/ 4 . It includes Hoyt 

(Vij = ^^ij = 0-5 ), and Nakagami-m {rjij = 1, fiij = m/2) as special cases. The analysis 
presented here can also be extended to Format 2. 

If a link is subjected to K-fi fading, the PDF of 7 *^ can be expressed as [l^ 



(3) 


where Kij > 0 denotes the ratio of the total power due to dominant components to the total 
power due to scattered waves. This fading includes Rice (/ij^ = 1 and Kij = K), and Nakagami- 
m {Kij — )■ 0 and fiij = m) as special cases. 

Let 7 i/i be the predetermined threshold for the sr link . If instantaneous SNR 7 ^^ of sr link 
falls below 7^/1, the link from node s to node d via node r is assumed to fail. On the other hand, 
if node r can successfully decode the messages from node s and forwards to node d, then the 
PDF of the end to end SNR 'ysrd can be given as [l29]l . (ISTIl 


flsM = ^r5(7) + (1 - 


(4) 


where Asr is the probability that outage occurs in sr link and 5(-) is Dirac delta function. The 
probability Agr in (H)) can be calculated as 


Agr 1 Prj^Ts?” ^ 'Jth} y 


(5) 


where Pr{-} is the probability operation and Fy^^{-) denotes the cumulative distribution of 
function of 7 ^^. 


III. Average Symbol Error Rate Analysis 


Mathematically, ASER P{e) of any modulation scheme for considered system can be computed 
as 03 


P(e) 



( 6 ) 
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where P{e\'y) is the eonditional SER of the modulation seheme in AWGN ehannels and f^Xl) 
is the the PDF of the total SNR 7 ^ of output of the MRC eombiner at node d. To obtain the 
elosed-form ASER expression, the moment generating funetion (MGF) of 7 ^ is required. The 7 * 
can be given as ||49l 


It = Isd + Isrd- (V) 

The MGF of 7 ^ is defined as Il49l 

poo 

= / exp(-z7)/.^„(7)d7, (8) 

Jo 

where u e {sd,rd, srd,t}, and ^ is the Eaplace variable. The of 74 can be obtained 

from the formula [|49ll 


~ -^7sd(^) -^7srd(^)- 

Using dH) in ([ 8 ]), the expression of can be obtained as ll43l 


(9) 


-^sr + (1 J^sr) ^ 

Finally, substituting (fT^ in dH), can be written as 


( 10 ) 


A^7j(^) Agr-M-y^^iz) {1 Asr)M'Yg^{z) . 


( 11 ) 


Using dll) in dS]), the (z) of 7-/1 faded ij link can be obtained as 

{2ijLij{hij + Hij) + 


MyY^) = 




Using dH) in dE]), the Ad.y. (z) of k-^ faded ij link can be obtained as Il43l 


MyY^) = 


f^ij (1 T ^ij ) 

f^ijY ^7) T 


^ij 


exp 




“ 1 “ ^ij) “ 1 “ ^'Jij 


( 12 ) 


(13) 


The scenario, where one hop of the link is subjected to rj-jd fading and the other hop is 
subjected to «:-/i fading, can be modeled accurately by mixed NFOS and EOS conditions, which 
occur in various applications including micro-/macrocellular and/or hybrid satellite/terrestrial 
communication systems mi, dn. We assume that all links experience independent and non 





identical but asymmetric fading. With two different fading types, there are six possible combina¬ 
tions for mixed fading scenarios as given in Table HI We have analyzed the mixed fading scenario 


TABLE I 

Fading parameters of dieferent scenarios. 


Scenario 

Links parameters 

sd 

sr 

rd 

1 

Vsd') fdsd 

^sri fdsr 

^rdi fdrd 

2 

^sdi l^sd 

VsE") 

'Qrd') l^rd 

3 

^sdi fdsd 

^sri fdsr 

^rd") fdrd 

4 

^sdi fdsd 

^sri fdsr 

^rd") fdrd 

5 

^sdi jdgd 

^sri fdsr 

^rdi fdrd 

6 

Vsd') jdsd 

^sri fdsr 

'^rd’} fdrd 


1 where 'jsd and 7 ^^ are modeled as rj-^ distribution but 'jrd is modeled as n-jj distribution. For 
rj-fi faded sr link can be given as [l47l 


A — 1 — V 

- -L //. 




ho 


‘^lisrhsr'Jth 

'Jsr 


(14) 


where Yiy{a, b) = —— ^ 

solution of this integral is ll47l 


Y„{a,b) = 1 - 


x^''q * I^^_i{ax^)dx is the Yacoub’s integral [HHl. A general 
i>]l + 2u] —(1 + a)P, —(1 — a)P) 


(1 - a2)-^6-4^r(l + 2z/) 


(15) 


where $ 2 ^^ is the confluent Lauricella function [|48l . In HTlI . another closed-form expression for 
the integer values of 2 /isr is also presented. 


A. Exact Analysis 

1) M-ary RQAM : The conditional SER performance of M = Mj x Mq-ary RQAM in 
AWGN channels is given as [|44|. Il50l 


F(e| 7 ) = 2 pQ (av^, 7 r/ 2 ) -f 2qQ {by/y,7i/2) - 2 pg {Q(a^, arccot( 6 /a)) 
+ Q ips/li arctan( 6 / a )}, 


(16) 
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where p = 1 — q = 1 — Mj and Mq are the number of in-phase and quadrature- 
phase eonstellation points, respeetively, a = ^ b = f3a and (] = dq/dj is the 

quadrature-to-in-phase deeision distanee ratio with dj and dq being the in-phase and quadrature 
decision distance, respectively, and Q{x,4>) is given as [321, 

Q{x,(j)) = -[ exp (--^-^)d9; a; > 0. (17) 

TT 7o ^ 2 Sin 6^/ 

Substituting (fT^ in ® and by algebraic manipulations ASER for RQAM scheme denoted as 
pRQAi^(e), can be expressed as 


FRQA“(e) = 2pX(a, 7r/2) -f 2gX(6,7r/2) — 2pg{X(a, arccot(6/a)) -f X(6, arctan(6/a)) }, 


(18) 


where the function X(-, •) is defined as 


X(a;,0) = / Qix^/^,(j))fy,{J)d-f 

Jo 

n oo 


exp 


2 

X 7 


2 sin^ 9 


fit{'l)dld9 


M 


It 


X 


2 sin^ 9 


d9. 


(19) 


Now, substituting (fTTI) in ([T9l) . X(x, 0) can be written as 


X(x,0) = AsrXi{x,(t)) -f (1 - A^^)X2 (x,0), 


( 20 ) 



(22) 
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The closed-form expressions for Xi(x, 0) and X 2 (a:, 0) are derived in Appendix. Thus, the ASER 
expression for given as 

FRQA“(e) = A5r{2pXi(a, 7r/2) + 2gXi ( 6 , 7r/2) — 2pgXi(a, arccot(6/a)) — 2pgXi ( 6 , arctan(6/a)) } 

+ (1 — A5r){2pX2(a, 7r/2) + 2gX2(6,7r/2) — 2pgX2(a, arccot(6/a)) — 2pgX2(6, arctan(6/a)) }. 

(23) 

For the special case of M-ary SQAmQ i.e. when Mj = Mq = \/M and 0 = 1, it can be 
shown that (1231) reduces to 

pSQAMj'g^ = 4A^^p{Xi(a,7r/2) - pXi (a, 7r/4)} + 4(1 - p{X 2 (a, 7r/2) - pX 2 (a, 7r/4)}, 

(24) 

where p = 1 — and a = ^ Foi" the special case of BPSK, i.e. when Mj = 2, Mq = 1 
and 0 = 0, it can be shown that (1231) reduces to 


FRRSR(e) = A,,Xi(x/2,7r/2) + (l-A,,)X2(v/2,7r/2). (25) 

2) M-ary XQAM : The conditional SER performance of M-ary XQAM in AWGN channels 
can be given as |[42l 

J /-1 

P{e\7) = glQ{ao^/y, 'k/ 2) -92'^ Q{ak\/l, arctan(A ;/{k + 1))) 

k=l 

V 

+ fi'2Q(ai\/7,7r/2) + 5(2 ^ arctan(A:/(A: - 1))) 

k=2 

v-l 

- fi'3Q(ao\/7,7r/4) - 25(2 ^ Q(aov^, arctan(l/(2A; + 1))), ( 26 ) 

k=l 

wher e M = 32,128, 512,..= 4 - ^2 = = 4 - ^ + f, z/ = ^, 

“0 = = V2kao, /c = 1, 2,..., Z/. 

Using a similar approach as followed to obtained (l23l) . a closed-form ASER expression of 


'Substitute y = x into ((42) and dmi to obtain Xi (a:,7r/4) and X 2 (a;,7r/4), respectively 
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XQAM can be obtained as given 


^xQAM(e) = (1 - A^^)|x2(ao,7r/2) + 5(2X2(01,7r/2) | + A^^|Xi(ao, 7r/2) 



k=2 


5 ( 3 Xi(ao, 7 r/ 4 ) - 2^2 (oo, arctan(l/(2fc + 1 


)))} + (1 - - 5(3X2(ao,vr/ 4 ) 


k=l 


V 


92 E^^ (ofc, arctan(fc/(fc — 1))) — 2^2 ^^X2(ao, arctan(l/(2fc + 1))) 


k=2 


k=l 



(27) 


k=l 

B. Asymptotic Analysis 

In asymptotic analysis, we assume that the average SNR values for three links are sufficiently 
large i.e., yij >> 1. With the high SNR assumption, the MGF of r]-p faded link in (fT2l) can be 
approximated as 



and the MGF of n-p, faded link in (fTSl) can be approximated as 



(29) 


Thus, Asr of sr link can be approximated as 



(30) 


where C is the inverse Laplace transform operator. For high SNR, we can also assume 

1. 
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1) M-ary RQAM: Using (l28l) . (I29l) and (1^ . the asymptotic ASER for RQAM, 
can be expressed as 

pRQAM(g) ^ yl^|2pX“(a,7r/2) + 2gX“(6,7r/2) - 2pgX“(a, arccot(6/a)) 

— 2pgX“(6, arctan(6/a)) } — 2pq'X“(a, arccot(6/a)) — 2pgX“(6, arctan(6/a)) 

+ 2pXr(a,7r/2) +2gX“(6,7r/2), (31) 


where 


xr 0 ) = 


TT 


M°° ( _ - _ 

V 2 sin^ 6 


de, 


(32) 


and 




M 


TT 


'ysd 


X 


2 sin^ 9 


M 


'^rd 




2sin^0 


dO. 


(33) 


The closed-form expressions for X“ (x, 0) and X“ (x, 0) are derived in Appendix. 

2) M-ary XQAM: Using an approach similar to that followed to obtain (l3TI) . a closed-form 
asymptotic ASER expression for XQAM can be given as 


v-l 

= A^|x“(ao,7r/2) - ^2 arctan(/c/(/c + 1))) 

k=l 

V 

+ 5f2X“(ai,7r/2) + 5(2 ^X“ (ofc, arctan(A;/(/c - 1))) 

k =2 

u-l 

- (ao, 7r/4) - 25(2 ^X“ (oq, arctan(l/ { 2 k + 1 ))) | 

k=l 

J/-1 

+ |x“(ao,7r/2) - 5(2 ^X“ (0^, arctan(/c/(fc + 1))) 

k=l 

u 

+ 5 ( 2 X 2 “(ai,7r/2) + g2^I^{ak,aTctan{k/{k - 1))) 

k =2 

1 /-I 

- 5 ( 3 X“(ao, 7 r/ 4 ) - 2 g 2 '^I^{ao,aTctan{l/{ 2 k + 1)))}. 

k=l 


( 34 ) 
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Fig. 1. ASER of 4 X 2-RQAM scheme with varying 77, k and fi. 


The asymptotic ASER expressions in (l3TI) and (1^ can be rewritten by substituting Pg = 
and Pr = {I — ^)P for 0 < .^ < 1, as 

where £i and £2 summarize the necessary elements in the asymptotic ASER expressions. The ^ 
is optimal power allocation factor. 


C. Diversity Order Analysis 

The diversity order, of dual-hop DE relaying system can be calculated as llSTTl 


C>df = lim 


P 

Nq- 


logpQAM^g^ 


(36) 


Using (l35l) and (|3^ . we can show that our considered scenario 1 (in Table H]) of this system 
achieves the maximum diversity order of C>df = ‘^dsd + min{2/4<jr, drd}- It can be observed that 
the fading parameters rfsdi Vsr and Krd have no impact on the diversity order. 
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Fig. 2. ASER of 32-XQAM scheme with varying rj, k and fi. 


D. Optimal Power Allocation Analysis 

Equal power allocation strategy for transmitters is not an optimal solution for power allocation. 
If partial channel state information is available at the transmitting nodes, an optimal power 
allocation can be performed to minimize the ASER at the destination node [|^ . The asymptotic 
ASER expressions in (|3TI) and (l34l) can be rewritten by substituting Pg = and = (1 “0-^ 
for 0 < .^ < 1, as 


pQAM(j) _ 


Q-^rsdQ-2fisrr^ n-‘^rsdc-2fi^an 
^‘‘sd ^‘‘sr ‘-'1 , “sd s ‘-'2 


+ 


(37) 


where Ci,C 2 > 0 are constants of the expressiorQ It can be shown that the second derivative of 
(1371) w.r.t ^ is always grater than or equal to 0, i.e., ^ 0- Hence (1371) is a convex 

function w.r.t. Equating the first derivative of (1371) w.r.t. ^ to zero, we get the relation 


-2Ci{psd + - {2psd + Prd)0- 


(38) 


^not given here as it is not required for optimization 
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Fig. 3. ASER of 4 X 2-RQAM scheme for dual-hop DF relaying system with varying and fixed rj = k = 1. 


It can be observed from (IMI) that the asymptotic optimum power allocation does not depend on 
the variance of link between s and d, it depends only on the variances of links between s and 
r and, between r and d. Moreover, we can also infer that the optimum transmitted power Pg at 
s is larger than 2jjisdP /+ j^rd) and less than P, while the optimum power Pr used at r is 
larger than 0 and less than firdP/i‘^fJ-sd + fJ-rd)- It means that we should always allocate more 
power at s and less power at r. 

IV. Numerical Results and Discussion 

Numerical examples concerning the ASER performance of QAM schemes are presented with 
computer simulations to verify the accuracy of analytical results. The Lauricella’s hypergeometric 
functions are numerically evaluated using their finite integral representation. In the numerical 
evaluation, we assume that the power is equally allocated to the source and the relay and the 
variance of each link is unity, i.e., ilgd = ^sr = ^rd = 1, unless stated otherwise. Additionally, 
the threshold 'jth is chosen according to jih = 2^^ — 1, with the spectrum efficiency R being set 
to 1 bit/s/Hz. 
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Fig. 4. ASER of 4 X 2-RQAM scheme for dual-hop DF relaying system with varying rj, k and fi. 


A. Direct Communication versus DF Cooperative Communication with RQAM and XQAM Schemes 

Figs.IH and [2] illustrate the ASER for 4 x 2-RQAM, and 32-XQAM sehemes, respeetively. We 
notiee that the numerieal results mateh with the eorresponding eomputer simulation results. These 
figures also show that the DF eooperative eommunieation ean always substantially improve the 
ASER performanee, relative to the direet eommunieation over the medium-to-high SNR region. 
For example, it ean be observed from Fig. [T] that for an ASER of 10“"^, P/Nq improvement that 
ean be aehieved in dual-hop DF relaying system {psd = f/sr = 1, Psd = Psr = 2, Krd = ‘2, Prd = 

2) over direet eommunieation (psd = 1, Psd = 2) is 3dB (approx.) 

B. Impact of Fading Parameters 

Fig. |3] shows the impaet of psd = Psr, and p^d for fixed rjsd = Vsr = x-rd = 1 on the ASER of 
4 X 2-RQAM seheme for dual-hop DF relaying system over mixed p-p and n-p fading ehannels. 

It ean be observed from this figure that the SNR gains are aehieved resulting from inereasing the 
fading parameters, psd = Psr, and p^d- For example, an ASER of 10“"^ with psd = Vsr = f^rd = ^ 
oeeurs at P/Nq ^ 34dB when psd = Psr = Prd = 0.5, P/Nq ^ 30dB when psd = Psr = 0.5, 
and Prd = 1, and, P/Nq ^ 26.5dB when psd = Psr = 1, and prd = 0.5. 
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Fig. 5. ASER of 4 X 2-RQAM scheme for dual-hop DF relaying system with equal and optimal allocation. 


Fig. |4] shows the impact of r], k and /i on the ASER of 4 x 2-RQAM scheme for dual-hop 
DF relaying system. Two important observations can be drawn from this figure. First, the SNR 
gains are resulting from increasing the fading parameters, k and /i, as expected. For example, an 
ASER of 10“'^ occurs at P/Nq ^ 24dB when rjsd = rjsr = 1, = ^J^sr = ^J^rd = 1, and Krd = 1, 

P/Nq 22dB when 7]sd = Vsr = 1, fJ-sd = fJ-sr = fJ-rd = 1, and Krd = 5, and P/Nq ^ 20dB when 
Vsd = Vsr = 1, ^J'sd = fJ-sr = fJ^rd = 2, and Krd = 1- Sccondly, there exists SNR penalty resulting 
with increase in the fading parameter, rj, as expected. For example, an ASER of 10“'^ occurs at 
P/Nq ^ 24dB when rjsd = Vsr = 1, V-sd = V-sr = l^rd = 1, and Krd = 1, and P/Nq 25dB when 

Vsd Vsr 5, V'sd V'sr V'rd !• 

C. Impact of Power Allocation 

Table |n] tabulates the fraction of the total power allocated to s for 4 x 2-RQAM as a function 
of Vtrd, Vsd = Vsr, Psd = Psr, K,rd and prd at P/Nq = 40dB. The values of the optimal power 
for any value of channel parameters of mixed p-p and K-p fading channels can be obtained by 
solving either (iTTl) or (IMl) using the software MATHEMATICA. Erom Eig. [51 we can see that 
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TABLE II 

Optimal Power Allocation Factor for 4 x 2-RQAM with /3 = 1, flsd = Osr = 1 and P/Nq = 40 dB. 


^rd 

Prd ~ 1 

f^sd f^sr 1 ? 

'^sd ~ ‘^sr ~ 1 

f^sd — l^sr 

^sd ^sr 


jdrd 

^rd 


1 

0.5 

0.01 

0.7760 

0.5 

1 

0.8000 

1 

0.6544 

4 

0.8000 

1 

0.01 

0.6706 

1 

1 

0.6682 

1 

0.6668 

4 

0.6679 

10 

0.5 

0.01 

0.9070 

0.5 

1 

0.8000 

1 

0.8228 

4 

0.8000 

1 

0.01 

0.6972 

1 

1 

0.6682 

1 

0.6682 

4 

0.6782 


when fading parameters and variances of link are highly unbalance, an optimal power allocation 
is useful for dual-hop DF relaying system. 

V. Conclusion 

Novel closed-form expressions are obtained for the ASER of general order RQAM and XQAM 
schemes with dual-hop DF relaying systems over mixed 17-/4 and /t-/i fading channels. The derived 
ASER expressions are used to analyze the performance of various QAM schemes with dual-hop 
DF relaying systems over mixed 17-/4 and n-jj, fading. We obtain the asymptotic ASER, which is 
useful to determine the system behavior at high SNRs in terms of the diversity order. Applying 
the asymptotic expression, optimal power allocation between the source node and the relay node 
under a total power constraint is analyzed. It is illustrated by numerical examples that the ASER 
performance vary according to the fading parameters of the communication links. All numerical 
results are validated through computer simulation results. 

Appendix 

Solution to the Integrals: Xi{x, 4>), l2{x, 0), I^{x, 0) and I^{x, 0) 

In the following subsections, we will derive the closed-form expressions for Xi{x, 0), X2{x, 0), 
X^{x,4>) and 0) in terms of multivariate Eauricella’s hypergeometric 
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functions which are given as [|45l . Il4^ 


and 


FP{a, 6i, 62, • • •, K] c; xi, X 2 ,,..., Xn) 


E 




m\,m2=0 


,,...,mn=0 * 2=1 


1 

B{a,c-a)Jo YY^=i{l - uXiY^ 
; Ixil < 1, \X 2 \ < 1 . . . \Xn\ < 1 


du 


(39) 


(a, 61,62, • • •, &n-i; c; xi, 0:2 ,,..., ai^) = 


B{a,c — a) Jq 


E ■ 

mi,m2=0„...,mn=0 


^{1—uY “ ^exp(Ma;n) 

\ 72 ivyi . 


du 


(c)Er 


n^.i+o) 


2=1 


where i?(a, 6) = is Beta funetion, mj = r(mi+l) and (a)n = r(a+ri)/r(a) 

is the Poehammer symbol for n > 0. The (•), and (■) funetions ean be easily and 
aeeurately evaluated by using its finite integral representation or its infinite series representation. 


A. Closed-form Expression for X\ (x, 7r/2) 

Substituting (fT^ into (1211) and followed by m = cos^ Q, the elosed-form expression for the 
integral Xi(x,7r/2) ean be given as ll43l 


X.P../2) ^ r(2;... + 0 5)>.,.,(xV2) 

2V?r(2^.i + i) 

^dsd(^hg(i T Hgfj \ 

^Psdifsd T Hgd') X 'ygd J 



^0.5, Pgdi Psd) ^Psd T 


^Esd{,dgd Hgd) 
‘^Esdi.hgd Hgd) A X '^gd 

(41) 


B. Closed-form Expression for Xi (x, arccot{y/x)) 

Substituting (fT^ into (l2TI) and followed by m = 1 — (y^/x^) tan^ 6, the elosed-form expression 
for the integral Xi(x, areeot(|//x)) ean be given as lf43l 


Xi (x, areeot(2//x)) = _ . „ E - Esd, Psd] ‘2pgd + 

27r(x^ + l/^)(2p^rf + 0.5) V 

^Psd{hsd Hgd) A X 'Jgd ^Pgd{hgd A Hgd) A X 'ygd \ 

4:pgd{hgd - Hgd) A (x2 A 2/^)%d’ 4:pgd{hgd A i^^d) A (x^ A y‘^)lsd) ' 


1.5; 


x^ 


x2 A ^2 ’ 


(42) 
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C. Closed-form Expression for X 2 (a:, 7r/2) 

Substituting, m and m into m followed by t = 
u = respectively, the closed-form expression for the integral X 2 (a:, 7r/2) can 

be given as 


X2(x, 7r/2) = 


X*!*! ( ‘2pg(i X prd X 0.5, 1, Psdi Psdi ‘^Psd X Prd X 1, " /i i ; \ ’ 

V X^lrd + ^Prd[l + Kd) 

‘^Prdif X krd'j'Jsd '^Psd^hgd Hsd'j'l/rd “^Prd^^ X krfj'ysd ‘^Psd^hgd Hgfj'Jrd 


exp( p^iiKrdfE f2pgd -\- prd X 0.5)x (2/ij.|^(l -|- rird)')^^'^{f^Pgd^J~h{p^rd/^ sd^^ 

2\/^ T(2pgd + Prd + l)(a:^7rd + 2/ird(l + 

2/^r(i(l X k^fj 


^ '^sd'^rd X 2,prd{f ~\~ krd)'^sd 
‘^Prd^i'dif X rird) 


^ '^sd'^rd X 2,prd{f krd)'^sd 


^‘^'1/rd X 2/ij.^(l -\- Krd) 


(43) 


D. Closed-form Expression for X 2 (x, arccot{y / x )) 

Substituting, m and m into m followed by t = and 

u = ^j t, respectively, the closed-form expression for the integral l 2 {x, arccot{y/x)) 
can be given as 


X 2 (^, arccot(^/a:)) = 


X$i j 2,fJjg(f -|- fJjpd H“ 0.5, 1, 0.5, flsd-i ^f^sd “1“ /^rd “1“ 1*5, 


6Xp( fdrd^rd^^'y/'Jrdi.^^f^sd i^rd “1“ 0.5) -|- 

27r {4.^sd^fh{^rdhsd))~‘^^^^{{x^ + |/^)7rd + 2/Xrd(l + 

2/^rd(l H“ ^rd) 


^ Trd “1“ 2/i;r^(l “1“ /^rd) 


^2 _|_ 2/X^rf(l + k 

rd) 

2/2-rd(l X krd'j'^sd ^Psdihsd Hgd)'^rd 


(x X y ')'yrd X 2^^|^(1 "1- Krd) X y )'ysd'Jrd X 2^,.£;(1 X krd^'Jsd 

‘^Prd(f ~k k^d)^gd ^Psd{kgd -\~ Hsd)'yrd “^Prd^^d^f X Krd) 


(x -\- y )'ygd '^rd X 2^j.£;(1 X krd)'ysd (a^ X y )'yrd X ‘2prdi^^ Krd) J 


(44) 


E. Closed-form Expression for Xi (x, aictanipj/z)) 

Substituting (fT^ into (l2T]) and followed by m = 1 — {z^ jy^) tan^ 6*, the closed-form expression 
for the integral Xi(x, aician{y/z)) can be given as ll43l 


Xi (x, arctan(?//^)) = 


y zM 


'fsd 


x'^iy'^+z'^) 

W 


:i(3) 


f ) f ) Psdi Psdi “^Psd X 1.5, 


y 


2Ti{y‘^ + z‘^){2pgd + D.h) ^ 

^Psd{hgd - Hgd)y‘^ + x'^y'^fgd 4:pgd{hgd + Hgd)y‘^ + x'^y'^fgd 


y^ + z^ 


^Pgd{hsd - Hgd)y^ + x 2(|/2 + z‘^)'jgd' Apgdihsd + Hgd)y^ + x^iy"^ + z‘^)'jgd 


( 45 ) 






















21 


F. Closed-form Expression for X 2 (x, arctan(|//z)) 

Substituting, m and m into m followed by t = 
u = rcspectivcly, the closed-form expression for the integral 

X 2 (x, arctan(|// 2 ;)) can be given as 


X 2 (x, arctan(?// 2 ;)) = 


X $1 ( 2,pgfi F Prd X 0.5, 1, 0.5, psdi Psd) ‘^Psd X Prd “fl.5, 2 \ 2^^ O fT ^ U ^ 2’ 

\ ^ yU “i~ ^ j'yrd “T -\- fv^djy 

^ y 'Jrd “1“ “h -|- kpd)y 'Jsd ^^sdi^hsd Hgd^y 'Jrd 


6Xp( yrd^rd)^yyj'^rdi^^sd f^rd 0.5) (2?/ “1“ 

27r {4y'^PsdVh{'yrd/%d))~‘^^‘’‘^{x‘^{y^ + z^)'yrd + 2y^prd{l + 

2prdi,^ X krd^y 


x2(^y2 z‘^)'yrd + 2prd{l + x‘^{y‘^ + z'^)'^sdlrd + 2prd{l + Kd)y‘^lsd 

2prd{l + krd)y'^fsd - ^Psd{hsd + Hgd)y'^frd 2pl^Krd{l + l^rd)y‘^ 


X2(^y2 + z‘^)'Jsd'lrd + 2prd{^ + krd)y‘^'lsd ’ + z‘^)'^rd + 2prd{l + l^rd)y‘^ 


(46) 


G. Closed-form Expression for If" (x,7r/2) 

Substituting (l28l) into (1^ and followed by m = cos^ 6, the closed-form expression for the 
integral Xf^(x, n/2) can be given as 


X“ (x, 7r/2, 'ysd) 


B{0.5,2pgd + 0.5) 

2n {4:PsdV^/{x^ Isd)) 


(47) 


El. Closed-form Expression for X|° (x,7r/2) 

Substituting, (l28l) and (l29l) into (1^ followed by m = cos^ 9, respectively, the closed-form 
expression for the integral X^(x, 7r/2) can be given as 

TOO I /rt\ ^ (0'6) 2pgd X prd X 0.5) exp( Prd^rd) [ /^rd(l X rird) 

^ - /T-/r 2- - ( -15—- 

2t: [Apgd\/hgd/{x^'^gd)) V x'^rd 



I. Closed-form Expression for Xf° (x, arccot{y/x)) 

Substituting (|2^ into (1^ and followed by m = 1 — {y‘^/x^) tan^ 9, the closed-form expression 
for the integral Xf^(x, arccot(?//x)) can be given as 


ir(..arcco.(,/.)) = 

x5(l,2/i,rf + 0.5)x(,')(^ 




xy 


27r(x^ -t- y"^) 
x^ 


1) 2psd X 1, 2pgd X 1.5, 


x2 _|_ y2 


(49) 
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J. Closed-form Expression for X|° (x, arccot{y/x)) 

Substituting, (1281) and (l29l) into (1^ followed by m = 1 — {y‘^/x‘^)tan^ 9, the closed-form 
expression for the integral X“(x, arccot(?//x)) can be given as 


^oo( tr / 'll f exp{-fj,rdi^rd)xyB{l,2psd + drd + 0.5) 

2 ix,arcco [y/x)) - + 27 r(x 2 + r/^) ((x^r/^)+/s:,,)))^- 

^ ^ ‘^dsd + Erd + 1; 2/isd + Prd + 1-5; ^ ) • (50) 


K. Closed-form Expression for Tf" (x, arctan(?// 2 ;)) 

Substituting (1281) into (1^ and followed by m = 1 — {y^ jz^) tan^ 9, the closed-form expression 
for the integral X“(x, axciaxi{y/z)) can be given as 


X“ (x, arctan(r// 2 ;)) = | 
X B (1, 2/i,d + 0.5) 


_4y^PsdVh^'\ zy 

x2[y2 z‘^)fsd) 27r(r/2 z"^) 

\ 1, ‘^Psd + 1; 2/isd + 1.5; —- 


(51) 


L. Closed-form Expression for X|° (x, arctan(?// 2 ;)) 

Substituting, (1281) and (l29l) into (1^ followed by m = 1 — {z^/y'^) tav? 9, the closed-form 
expression for the integral X^(x, aician{y/z)) can be given as 


X“ (x, arctan(r// 2 ;)) = 


X Fj^ ( 1, 2pgii F prd X 1, 2pg(i X prd F 1.5, 


exp{-prdi^rd)z yB (1, 2psd + Prd + 0.5) 


^y’^psdV^ _ 

x2(^2 ^ Z^)%dJ 27r(y^ + Z^)(x^(y^ + Z^)frd/(y^Prd(l + 

2 


y 


y^ + 


(52) 
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